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Abstract 

We study the Zee model of neutrino mass in the framework of i?-parity 
violating supersymmetry. Within the matter contents of the minimal super- 
symmetric standard model, any one of the three right-handed sleptons could 
be a suitable candidate for the charged-singlet scalar of the Zee model, and 
one of the Higgs doublets provides the extra necessary vacuum expectation 
value. A combination of one bilinear and two trilinear 72-parity-violating cou- 
plings then completes the model. In this framework, we also discuss other 
various contributions to neutrino masses and derive the conditions for the 
dominance of the contribution from the Zee model, and hence maintain the 
successfully Zee mass texture. However, this model within the minimal su- 
persymmetric standard model is shown to be only marginally feasible. More 
general versions of supersymmetrization of the Zee model are also discussed. 
A particularly interesting example that has extra Higgs doublets while the 
slepton, especially the selectron, keeps the role of the Zee scalar is illustrated. 
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I. INTRODUCTION 



Within the standard model (SM), the V — A nature of the weak interaction dictates a zero 
mass for all three families of neutrinos. The discovery of neutrino mass(es) or oscillation(s) 
will certainly push for new physics. Evidence for neutrino oscillations has been collected 
in a number of solar neutrino and atmospheric neutrino experiments. The most impressive 
results were the recent neutrino deficit and the asymmetric zenith-angle distribution 
observed by the Super-Kamiokande Collaboration 

The atmospheric neutrino deficit and zenith-angle distribution can be explained by the 
Vfj, — v T or z/ M — v s oscillations [y s is a sterile neutrino that has negligible coupling to W or 
Z boson and the former has a slightly better fit). The oscillation parameters with — > v T 
at 90% C.L. are [| 

Ami - (2 - 6) x 10- 3 eV 2 , sin 2 2fl atm > 0.85 . 

On the other hand, the solar neutrino deficit admits more than one solution. With v e — > v T 
the solutions at 95% C.L. are @ 

vacuum oscillation: Am 2 ol ~ (5 - 8) x 10 -11 eV 2 , sin 2 20 so i ~ 0.6 - 1.0 , 
small angle MSW: Am 2 ol ~ (4 - 9) x 10" 6 eV 2 , sin 2 26> so i ~ (3.5 - 13) x 10~ 3 , 
large angle MSW: Am 2 ol ~ (8 - 30) x 10" 6 eV 2 , sin 2 29 sol ~ 0.4 - 0.8 . 

The above two neutrino-mass differences can be accommodated by the three species of 
neutrinos that we know from the SM. There is also another indication for neutrino oscillation 
from the accelerator experiment at the Liquid Scintillation Neutrino Detector (LSND) [BJ, 
which requires an oscillation of into another neutrino with 

Am 2 SND ~ 0.2 - 2 eV 2 , sin 2 2# LSND ^ 0.003 - 0.03 . 

To as well accommodate this data it requires an additional species of neutrino beyond the 
usual neutrinos. Nevertheless, further evidence from the next round of neutrino experiments 
is required to confirm the neutrino mass and oscillation. As neutrino is favored to be massive 
it is desirable to understand the generation of neutrino masses from physics beyond the SM, 
especially, to see if the new physics can give a neutrino mass pattern that can explain the 
atmospheric and solar neutrino data, and perhaps the LSND as well. 

An economical way to generate small neutrino masses with a phenomenologically favor- 
able texture is given by the Zee model |§-§| , which generates masses via one- loop diagrams. 
The model consists of a charged gauge singlet scalar h~, the Zee scalar, which couples to 
lepton doublets ip Lj via the interaction 

f (VfiCVfi) W h~ , (1) 
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where a, (3 are the SU(2) indices, i,j are the generation indices, C is the charge-conjugation 
matrix, and / y ' are Yukawa couplings antisymmetric in i and j. Another ingredient of the 
Zee model is an extra Higgs doublet (in addition to the one that gives masses to charged 
leptons) that develops a vacuum expectation value (VEV) and thus provides mass mixing 
between the charged Higgs boson and the Zee scalar boson. The corresponding coupling, 
together with the / y 's, enforces lepton number violation. The one-loop mechanism for the 
Zee model can be found in Fig. [I]. 

A recent analysis by Frampton and Glashow || (see also Ref. ||) showed that the Zee 
mass matrix of the following texture 
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where e is small compared with m efl and m er , is able to provide a compatible mass pattern 
that explains the atmospheric and solar neutrino data. The generic Zee model guarantees 
the vanishing of the diagonal elements, while the suppression of the m^ T entry, here denoted 
by the small parameter e, has to be otherwise enforced. Moreover, m efJi ~ m eT is required 
to give the maximal mixing solution for the atmospheric neutrinos. We shall describe the 
features in more detail in the next section. 

So far the Zee model is not embedded into any grand unified theories or supersymmetric 
models. Here we analyze the embedding of the Zee model into the minimal supersymmetric 
standard model (MSSM) with minimal extensions, namely, the .R-parity violation. The right- 
handed sleptons in SUSY have the right quantum number to play the role of the charged Zee 
scalar. The .R-parity- violating A-type couplings (XLLE) could provide the terms in Eq.(^). 
It is also easy to see that the i?-parity-violating bilinear //-type couplings (fiiLH 2 ) would 
allow the second Higgs doublet H 2 in SUSY to be the second ingredient of the Zee model. So 
far so good. However, the SUSY framework dictates extra contributions to neutrino masses, 
which deviate from the texture of the Zee mass matrix of Eq. (0). The major objective 
of this paper is to address the feasibility of the embedding and to determine under what 
conditions could one make a supersymmetric Zee model within the -R-parity- violating SUSY 
framework while retaining the successful flavor of the former. We will also discuss briefly 
more generic versions of supersymmetric Zee model. There is also a study of Zee mass 
matrix within the framework of gauge-mediated SUSY breaking with the messenger field as 
the Zee singlet [|7|. 

In .R-parity- violating SUSY, there are three other sources for neutrino masses, in addition 
to the Zee model contribution. They are (i) the tree-level mixing with the higgsinos and 
gauginos, (ii) the one-loop diagram that involves the usual mass mixing between the left- 
handed and right-handed sleptons proportional to m? (Af — /itan/?), and (iii) the one-loop 
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diagram that again involves the mixing between the left-handed and right-handed sleptons 
but this time via the A and fa couplings. [] The first two contributions have been considered 
extensively in literature ||, but the last one is identified here for the first time. Also, we 
are the first one to identify the Zee model contribution to the neutrino mass in the SUSY 
framework. Furthermore, we will obtain the conditions for the Zee model contribution to 
dominate over the contributions in (i) and (ii). The contribution in (iii) can actually preserve 
the texture of the mass matrix of Eq. 

There are complications in choosing a flavor basis when R parity is broken. Actually, 
the form and structure of the lepton mass matrices under the coexistence of bilinear and 
trilinear i?-parity- violating couplings are basis dependent. One has to be particularly careful 



with a consistent choice of flavor basis. Here we adopt the single- VEV parametrization flO 



that provides an efficient framework for our study. The most important point to note here 
is that this parametrization implies a choice of flavor basis under which all three sneutrinos 
have no VEV, without any input assumptions. All i?-parity-violating couplings introduced 
below are to be interpreted under this basis choice. 

The organization of the paper is as follows. In the next section, we describe the texture 
of the Zee mass matrix in Eq. @. In Sec. Ill, we calculate entries in the 3x3 neutrino 
mass matrix from all the sources listed above. In Sec. IV, we derive the conditions for 
the contributions from the Zee model and from (iii) above to be dominant. Section V is 
devoted to discussions on more general versions of supersymmetrization of the Zee model. 
We conclude in Sec. VI. 



II. ZEE MASS MATRIX 



Here we briefly describe the basic features of the Zee mass matrix, as given in Eq.(^). 
We first take e = 0. The matrix can be diagonalized by the following transformation 
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(3) 



1 There is one other type of contribution from a gaugino-sneutrino loop with neutrino-antineutrino 
mass splitting. The latter could be a result of i?-parity violation. This was discussed in Ref. ||. 
This contribution depends on the -R-parity violating B terms and basically enters in the same 
entries in the neutrino mass matrix as the corresponding m terms, and thus supplementing the 
latter. We will neglect this type of contribution in this paper. 



4 



with the eigenvalues m, — m, for v L ^v h ^v h ^ respectively, and m = ^m 2 eil + m 2 er . Hence, 
the two massive states form a Dirac pair. The atmospheric mass-squared difference Am 2 tm ~ 
3 x lCT 3 eV 2 , is to be identified with m 2 = m 2 + m 2 T . The transition probabilities for v Lfi 
are 



D = 



m P „m 



P, Jr ->„. = 4 °" *' sin 



efi"'er 1 ■21 \ e/u 1 er 



mi, + mL)L 



If m etl ~ m er , then sin 2 2# atm ~ 1. This mixing angle is exactly what is required in the 
atmospheric neutrino data. The neutrino mass matrix texture with e = can be called 
the zeroth order Zee texture. It is the first thing to aim at in our supersymmetric model 
discussions in the next section. 

If we choose a nonzero e, but keep e <C rn eil ^ eT . Then after diagonalizing the matrix we 
have the following eigenvalues 



m u2 = -dm? + m T + e^— , 

V m 2 efl + m 2 eT 

m„ 3 = — 2e- 



The mass-square difference between vr? ul and m 2 2 can be fitted to the solar neutrino mass. 
For instance, one can take the large angle Mikheyev-Smirnov-Wolfenstein (MSW) solution 
and requires 



giving 



At me ^ meT = Am s 2 ol ~ 2 x 1(T 5 eV 2 
Imlp + m 2 er 



5 x 10 -3 



where we have used m e[l ~ m eT . 

We will see below that in our supersymmetric model the couplings that are required to 
generate the zeroth order Zee texture also give rise to other contributions, which have to be 
kept subdominating in order to maintain the texture and hence the favor of the Zee model. 
Even though these extra contributions might not be identified as the same entries as the e 
parameter of Eq. @, i.e., appear in the diagonal entries instead, they could still play the 
same role as to give a phenomenologically viable first order result for a modified Zee matrix. 
Hence, we will not commit ourselves to the first order Zee matrix as given in Eq.(^), but only 
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to its zeroth order form, namely with e = 0. The first order perturbation is then allowed 
to come in through any matrix entry. It will split the mass square degeneracy of the Dirac 
pair similar to the e case above. For example, if the first order perturbation is given by a e d 
appearing at the m ee entry, the resulting mass eigenvalues are modified to 

m vl = ^m 2 efl + m 2 eT + ^ , 
m u3 = . 

The mass-square difference between m 2 1 and m 2 2 can then be fitted to the solar neutrino 
data and we obtain e d /m eM ~5x 1CT 3 , the same as above. If, on the other hand, e d appears 
at the or m TT entry, the solar neutrino data can still be fitted and t d /m efJl ~ 1 x 1CT 2 
is required. Once the zeroth order Zee texture for the atmospheric neutrino is satisfied, it 
is straightforward to further impose the above condition for the solar neutrino. 

III. NEUTRINO MASS MATRIX 

First consider the superpotential as given by 

W = eJ^CRHfU? + Y?Q?HfDj + YZLfHfE? + pH?H? + \ jk Lt^E c k + ^LfH 2 ^ , 

(4) 

where e 12 = — e 21 = —1, i,j = 1, 2, 3 are the generation indices. H 1 = (h®, h~), H 2 = (hf, h° 2 ). 
In the above equation, Q, L, U c , D c , E c , H u and H 2 denote the quark doublet, lepton doublet, 
up-quark singlet, down-quark singlet, lepton singlet, and the two Higgs doublet superfields. 
Here we allow only the i?-parity violation through the terms LLE C and LH 2 with coefficients 
Ajjfc (antisymmetric in and respectively. The other i?-parity-violating couplings are 
dropped as they are certainly beyond the minimal framework needed for embedding the Zee 
model. The soft SUSY breaking terms that are relevant to our study are 

(Y^ijLfHfE? + (\A x ) ljk LtL^E c k + pBHfH? . 

Actually, the (XA X ) terms do not contribute because our choice of basis eliminates the 
VEV's for Lj's. This simplifies the analysis without lose of generality. We adopt the single- 
VEV parametrization, which uses the L t basis such that the charged-lepton Yukawa matrix 
Y E is diagonal. The whole (Y E A E ) term will be taken as predominantly diagonal, namely, 
(Y E A E )ij ks Y E A E 8ij (no sum). This is just the common practice of suppressing off-diagonal 
A terms, favored by flavor-changing neutral-current constraints. 



6 



The tree-level mixing among the higgsinos, gauginos, and neutrinos gives rise to a 7 x 7 
neutral fermion mass matrix M m- 
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whose basis is {—iB 1 —iW,h^ ) h° 1) i , La) v l 



\ 



(5) 



i ,i4 r ). Each of the charged-lepton states deviates 
from its physical state as a result of its mixing with higgsino-gaugino through the corre- 



sponding \ii term p0| . However, we are interested only in a region of the parameter space 



where the concerned deviations are negligible, as also discussed in Ref. [11]. Hence, we are 
effectively in the basis of the physical charged-lepton states, as indicated. In the above 
7x7 matrix, the whole lower-right 3x3 block (m") is zero at tree level. They are induced 
via one-loop contributions to be discussed below. One-loop contributions to the other zero 
entries are neglected. We can write the mass matrix in the form of block submatrices: 



Mm 



M 









(6) 



where M is the upper-left 4x4 neutralino mass matrix, £ is the 3x4 block, and m° v is the 
lower-right 3x3 neutrino block in the 7x7 matrix. The resulting neutrino mass matrix 
after block diagonalization is given by 



<M- % e + (ml) 



(?) 



The first term here corresponds to tree level contributions, which are, however, see-saw 
suppressed. 

Before going into our best scenario analysis, we will sketch how the couplings, A^'s and 
/ii's, lead to the neutrino mass terms. While some of them have been studied in litera- 
ture, others are identified here for the first time. We do this from the perspectives of the 
supersymmetric Zee model, but the results are quite general. 

Our minimalistic strategy says that a A^ or a fii should be taken as zero unless it is 
needed for the Zee mechanism to generate the neutrino mass terms m eA1 and m eT . Readers 
who find the extensive use of unspecified indices in the following discussions difficult to 
follow are suggested to match them with the results for the explicit examples that we will 
list below. We identify the following four neutrino-mass generation mechanisms. 
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(i) Zee mechanism. We show in Fig. [T], the two Zee diagrams for the one-loop neutrino 
mass terms. The right-handed slepton £ Rk is identified as the charged- singlet boson of the 
Zee model, and its coupling to lepton fields has the correct antisymmetric generation indices: 
see Eq. (JJ). To complete the diagram the charged Higgs boson h\ from the Higgs doublet 
H 1 is on the other side of the loop and a iR k -h[ mixing is needed at the top of the loop. Such 
a mixing is provided by a F term of Lk- n^vrkjij^ (h®)/(h®), where h® takes on its VEV, 
for a nonzero /ifc. Thus, the neutrino mass term (m°)y has a 

^kTTTt k \jk{mf. -vrQ (8) 

dependence, where r/^.'s are the charged lepton masses. 

(ii) LR slepton mass mixing. Another, well-studied, type of contributions comes from 
the one-loop diagram with two A-coupling vertices and the usual (A E — // tan /3)-type LR 
slepton mixing. Neglecting the off-diagonal entries in A E , the contribution to (m")y with 
the pair A^ and Xjki is proportional to 

[ (A E k - n tan/3) + (1 - 8 H ) (A E - //tan/3) ] m ek m ei \ Uk \ jkl . (9) 

Only Ajjfc's with all distinct indices (e.g. X m ) fail to give contributions of this kind on its 
own. A nonzero A^ contributes to the diagonal (m°)jj. An illustration for the term is given 
in Fig. |2|. With any two nonzero A^'s, this kind of contributions to the (mj) entries, in 
particular the diagonal ones, cannot be avoided. 

(iii) LR slepton mass mixing via R-parity violating couplings. This contribution is iden- 
tified here for the first time. While the contributions to generation mixing in the usual 
(A E — yutan/3)-type LR slepton mixing via the off-diagonal entries in A E are expected to be 
small, there is another independent source of generation mixing in the LR slepton-mass mix- 
ing, which may not follow the rule. The latter comes from a F term of Lf. iMKjk^^R (h®), 
where h® takes on the VEV. This is similar to the £a,-h~ mixing in the Zee model, except that 
this time we have a A-type coupling instead of the i?-parity-conserving Yukawa coupling. 
This newly identified source of mixing results in constraints on the ^iXijk products, which 
is an interesting subject of lepton-slepton phenomenology studies. 

With a specific choice of a set of nonzero /ij's and A's, this type of mixing gives rise to 
the off-diagonal (mj)q terms only and, therefore, of particular interest to our perspectives 
of Zee model. Taking the pair A^ and Xjhi for the fermion vertices and a F term of L g 
providing a coupling for the scalar vertex in the presence of a \i g and a X g hk (see Fig. |[), a 
{wnfyij term is generated and proportional to 

H>g m i l ^ghkXukXjhl ■ (10) 

The proliferation of indices here is certainly difficult to keep track of. When we allow only 
a single nonzero A at a time, the only contribution comes from Ayj but not from those with 
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distinct indices. Suppose we have nonzero Xyj and fa, they then give a contribution to the 
off-diagonal {m° y )ij with a 

/'./"'.. 

dependence, which is obtained from expression ( |XD|) through the substitution h = i and 
g = k = I = j. It is easy to see that for a minimal set of nonzero fa and A^ required 
to generate the zeroth order Zee texture, this minimal set also contributes to the same 
neutrino mass terms via the new mechanism identified here. Hence, they are desirable from 
the perspectives of keeping the Zee mass matrix texture. 

(iv) Tree-level mixing. Through gaugino-higgsino mixings, nonzero fa's give tree-level 
see-saw type contributions to (m v )ij proportional to fa fa, i.e., through the first term in 
Eq.(^) instead of the second. With the contribution put in explicitly, Eq.(|7|) then gives 



v 2 cos 2 /? (g 2 M 1 + g' 2 M 2 ) 

~ ~ 2/i \2faMJM 2 - v 2 sin/? cos/3 {g 2 M 1 + g' 2 M 2 )} 



(III,,),, = - - — — — — : — r— — — r fa fa + ('/'"),, • (11) 



A diagonal (m v )kk term is always present for a nonzero fa as needed in the Zee mechanism. 
This contribution has no charged lepton mass dependence. To eliminate these tree-level 
terms requires either very stringent constraints on the parameter space or extra Higgs su- 
perfields beyond the MSSM spectrum. We will see that this is a major difficulty of the 
present MSSM formulation of supersymmetric Zee model. 

From the above discussions, we conclude that a minimal set of i?-parity violating cou- 
plings needed to give the zeroth order Zee matrix is the following : 

{ A 12 fc , A 13 fc , fa } • 

As at least one of the two A's has the form A^ (= —Xkik), all types of contributions that 
have been discussed above are there. We want to make the contribution from the Zee 
mechanism dominate over other contributions, or at least to make the diagonal mass entries 
to (m u ) subdominant. This necessarily requires subdomination of the contributions from 
the tree-level see-saw mechanism and from the (A E — fi tan/3) -type LR slepton mixing. So, 
it is the Zee mechanism and the newly identified mechanism, which involve the interplay 
between the bilinear fa and trilinear A^ i?-parity-violating couplings, that are required to 
be the dominating ones. 

We will discuss below two illustrative scenarios (1) X m , A m , and \i x and (2) A E3 , A B 3, 
and After the 7x7 matrix M.m is block diagonalized, the resulting 3x3 neutrino mass 
matrix (m„) of Eqs.(|7|) or ([11]) is obtained in each of these scenarios. 

Scenario 1: \ m , \ m , and fa. The resulting neutrino mass matrix for scenario 1 is given 

by 
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( C-tnl C 2 m e rn 2 l ^ 1 \j A +C A rn e ^ 1 \ 1A C 2 m e m 2 T fi 1 \ m + C 3 m e fiA m ^ 
(7 4 m 2 A^ 2C 4 m 2 A^i 



V 



C 4 m 2 



(12) 



which is symmetric and we only write down the upper triangle. The Cj's are given by 

v 2 cos 2 (3 (g 2 M 1 + g' 2 M 2 ) 



c 2 
c, 
c 4 



2/i [2fxM 1 M 2 - v 2 sin/3 cos/3 (g 2 M 1 + g' 2 M 2 )] ' 
-]_V2 tan/3 

16tt 2 u cos/5 Jli %' iWe fi J ' 



1 f sin/3 
16tt 2 
1 



AL/^MfJ + A^/^MfJ 



16tt 2 



(^-/x tan/3) f(M~ L , M 2 ) , 



(13) 



where 



x-y 



log 



Scenario 2: \ n3 , A B3; and /i 3 . For scenario 2, the neutrino mass matrix is given by 



< C' 4 m 2 A| 3 C' 2 m T m 2 /^A^ + C 5 m T /'. ; A,,..A-',. C2 ^ 3 A B3 + C B m T MsA| 3 N 


V u 2 / 



(14) 



where 



C' = 



16vr 2 



tan/3) f(M 2 L ,M 2 R ), 



-1 v 7 ^ tan/3 



a = ^A^^/(M 2 M? R ), 



a 



16n 2 v cos/? 
1 t> sin/3 



167T 2 



f(M 2 L ,M 2 R ) 



(15) 



In the above, we have neglected terms suppressed by m e /m M or m e /m T . There is also another 
scenario, with {A K2 , A 132 ,/i 2 }, which is very similar to this scenario 2. 



IV. CONDITIONS FOR MAINTAINING THE ZEE MASS TEXTURE 

In order to maintain the zeroth order Zee texture as discussed in Sec. II, we need m e[ii 
and m eT to dominate over the other entries. Moreover, we need m eM ~ m eT ~ yAM 2 tm (~ 
5x 10~ n GeV). Here we give an estimate of the required conditions on the model parameters, 
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for a chosen minimal set of i?-parity-violating couplings {A 12 fc, Xi 3 k, Hk] (with a specific k)}. 
Since we are interested only in the absolute value of each term and so we will drop negative 
signs wherever feasible. We will look at each matrix entry in Eqs. flT2|) and ([H]) carefully. 

Scenario 1. Requiring the tree-level gaugino-higgsino mixing contribution to m ee in Eq. 
(|12|) to be negligible compared to m efJi , it gives 

[i\ cos 2 /3 < jj 2 M 1 (1 x 1(T 14 GeV" 1 ) , (16) 

in which we have assumed M x « 0.5M 2 . The condition is rather stringent that requires 
either very small fa, fa cos/3 <C 10~ 4 GeV at M^/i ~ 0(100) GeV, or particularly large 



gaugino mass(es). As pointed out in Ref. [TO], the dependence on tan/5 is very important 



here. The cos/3 goes from order one to ~ 0.02 in the domain of large tan/3. 

Requiring the (A1 — fi tan/3) LR slepton mixing contribution to be much smaller than 



m eil , we have 



„ „ max(M| , Mf ) , „ , , 

AL, A 2 m « £L °y (3 x 10- 2 GeV" 1 ) , (17) 



— /i tan/3) 

where we have used 



f(M*Ml) ~max(M^,M4). (18) 



The constraint in Eq. (|17D is obviously very weak. In fact, it can certainly be neglected, 
especially when other phenomenological constraints |6|;[l2j on A's [as effective Zee couplings 



fij of Eq. (|1|)] are taken into consideration. The relevant constraint here is given as 

A? 



m < 10- 8 GeV" 2 , (19) 



M- 



from the tree- level Zee-scalar mediated /i decay jl^]. The upper bound on A^ is hence no 
better than 0.01 for Mg R at 100 GeV. The corresponding constraint on A^ from r decay is 
definitely weaker, which has no relevance here as we will see below. Hence, the suppression 
needed for m m , m^ T , and m TT of Eq. ( |T2"D is easy to obtain. 

The remaining question is if one can still generate the right (order of) m efJi and m er 
when fa, X m and A^ satisfy the above constraints. Let us first look at the m e[li entry. From 
Eq . (|12"1) , m efJL has two contributions. The first one (the one with a C 2 dependence) is from 
the authentic Zee mechanism. For this contribution to give the right m eil value, it requires 

m *> ~ ^ (2 x 10- 10 GeV 2 ) ~ (5 x lO" 11 GeV) , (20) 



cos 2 /3 max(M 2 _,Mf 



2 



or 



(fa cos/3) \ m ~ cos 3 /3 max(M 2 _, MfJ (0.25 GeV" 1 ) . (21) 



11 



The right-hand side above cannot be much smaller than 0.01 GeV, even at the more favorable 
case of very large tan/3. Given the above constraints in Eqs. ([T6| ) and (|19|) , this is obviously 
unrealistic. Even though the corresponding contribution (with a C 2 dependence) to m eT has 
a m 2 /m 2 enhancement and depends on A^ instead of X m , it is not much better than m e/Ji . 

For the second contribution (the one with a C 3 dependence), to give the right m e/1 value, 
it requires 



\ 3 \ 2 \ 



max(M? L ,M|J max(M3, Mf^, } 



(5 x lO^GeV 2 ) ~ (5 x 10~ n GeV) . 

(22) 



A naive comparison with Eq. ( ^0|) above illustrates one important fact. Assuming a common 
scale for the scalar masses, the A coupling(s) only have to be larger than 10~ 3 for this second 
contribution to be larger than the first one. From Eqs. (O) and (|T9|), such A's are easily 
admissible. With X m « A^ the corresponding contribution to m er has the same form as 
m eA1 , with the interchange of A^ with X m , and thus is of a similar value. The condition in 
Eq. (|22"D then becomes 



(/ij cos/3) \ Ul 



where we have taken Mg R to be the dominating mass among the scalars. The latter choice 
corresponds to the optimal case because smaller scalar masses help reducing the size of the 
Aiii needed, while on the other hand, in Eq. fll9|) larger Mg R relaxes the constraint on Af 4l . 
With Eqs. (|IED and (p!9| ) taken into consideration, the result ends up actually no better than 
the best (large tan/3) case of Eq.(pO|) above. 

Scenario 2. Here we follow our above analysis for this more interesting scenario. Requir- 
ing the tree-level gaugino-higgsino mixing contribution to be well below m e/1 gives 

Hi cos 2 /3 < /i 2 M x (1 x 10- 14 GeV- 1 ) . (24) 

This is basically the same as in scenario 1, though it corresponds to m TT instead. 

For the (Af; — fi tan/3) LR slepton mixing contribution to be much smaller than m eA1 , we 
have 

Max(M? ,M? ) , 

A °-« i4--^ '" xll M ) ' (25) 

This corresponds to m ee . It tells us that A B3 can hardly be much larger than 10~ 3 . On 
the other hand, A^ is constrained differently because it does not contribute to this type of 
neutrino mass term. The constraint that corresponds to Eq . (|l~9D , however, becomes 

^<10- 8 GeV" 2 , (26) 
M? R ~ 
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which tells us that A^ can be as large as order of 0.01 for scalar masses of order of 0(100) 
GeV. 

Again both m eM and m er have two terms. Let us look at m efJi first. For the first term in 
m eA1 (the one with a C 2 dependence) in Eq. ( |i4|) to give the required value of atmospheric 
neutrino mass, we need 

meu ~ TTT-o » ,9 x (7 x lO^GeV 2 ) ~ (5x lO" 11 GeV) (27) 

M cos 2 /3 max M 2 M| v ; v ; K J 

v h 1 R 

or 

(jjl 3 cos/5) A]23 ~ cos 3 /? max(M h 2 -, M? R ) (7 x 10~ 5 GeV" 1 ) . (28) 

This result looks relatively promising. If we take cos/3 = 0.02, all the involved scalar masses 
at 100 GeV and A E3 at the corresponding limiting 0.01 value, /x 3 cos/3 has to be at 5.6 x 
10~ 4 GeV to fit the requirement. This means pushing for larger M 1 (and M 2 ) and fi values 
but may not be ruled out. 

What about the corresponding first term in m eT entry? The term has a A B3 dependence 
in the place of A^ with an extra enhancement of m 2 /m 2 , in comparison to m efl . That is to 
say, requiring m eA1 « m er gives, in this case, 

2 

Al33 ~ n A]23 • (29) 

This gives a small A B3 easily satisfying Eq. (P5p. The small A B3 also suppresses the second 
terms in both m eM and m er , the C 5 dependent terms in Eq. flU]) . Note that the above 
equation represents a kind of fine-tuned relation between the two couplings A B3 and X n3 . 
More precisely, the value of Aj3 3 m 2 has to be within a factor of 1.8 of that of A^ m 2 in order 
to fit sin 2 26 l atm . This feature is inherited directly from the original Zee model, as discussed 
in Ref. 0. Nevertheless, it is difficult to motivate this relation from a theoretical point of 
view. Phenomenologically, the relation implies that A B3 is two orders of magnitude smaller 
than Aj2 3 , which indicates a strongly inverted hierarchy against the familiar flavor structure 
among quarks and charged leptons. Since the current experimental bounds from the rare 
processes, such as //, r — > ej, showed the usual hierarchical trend down the families, the 
relation in Eq. (p9| ) says that once the constraints on the A^ are satisfied, A B3 should be 
automatically safe. This justifies our above statement that the r-decay constraint analogous 
to Eqs. ( |T9"D and ( |26| ) have no relevancy here. Conversely, if A B3 contributions to some rare 
processes are identified in the near future, it would spell trouble for the SUSY Zee model 
discussed here. 

Finally, we comment on whether it is feasible to have an alternative situation in which the 
second (C s dependent) terms in m eA1 and m eT dominate over the first (C 2 dependent) terms. 
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The comparison between these two types of contributions is similar to that of scenario 1, 
as can be easily seen by comparing terms in Eqs. ([12]) and (0). As in scenario 1, we need 
to push A B 3 to the order of 0.01. This at the same time requires either a particularly large 
Mf L or some fine-tuned cancellation between A® and [i tan/3 in order to fulfill the condition 
in Eq. fl25|) . Thus, it is unlikely to have the second terms of m eil and m eT dominant over 
the first terms. 

To produce the neutrino mass matrix beyond the zeroth order Zee texture, the subdom- 
inating first-order contributions are required to be substantially smaller in order to fit the 
solar neutrino data. Here, it is obvious that it is difficult to further suppress the tree level 
gaugino-higgsino mixing contribution to m TT , which makes it even more difficult to get the 
scenario to work. Explicitly, the requirement for the solar neutrino is 

l4 cos 2 /? ~ ^M, (1 x 10~ 16 GeV^ 1 ) , (30) 

following directly from the result given in Sec. II [cf. Eq. (f24|)]. 



V. MORE GENERAL VERSIONS OF SUPERSYMMETRIC ZEE MODEL 

We have discussed in detail the minimalistic embedding of the Zee model into the minimal 
super symmetric standard model. The conditions for maintaining the Zee neutrino mass 
matrix texture is extremely stringent, if not impossible. Here we discuss some more general 
versions of supersymmetrization of the Zee model. 

As mentioned in the Introduction, an easy way to complete the Zee diagram without 
the /Xj-type, bilinear _R-parity-violating, couplings is to introduce an additional pair of Higgs 
doublet superfields. Denoting them by H 3 and H 4 , bearing the same quantum numbers as 
H 1 and H 2 , respectively, i?-parity- violating terms of the form 

can be introduced. With a trivial extension of notations (in Fig. |l]with replaced by h^), 
we obtain a Zee diagram contribution to (m, v )ij through as follows : 

Tei ft K - O x «* ^ K M l v M l k ) • (3i) 

Here the slepton l Rk keeps the role of the Zee scalar. We have neglected the F term obtainable 
in the existence of bilinear \i terms between H 3 and H 2 or H 4 . At least one of them has to 
be there. When no LiH 3 Ej type i?-parity-conserving Yukawa couplings are allowed, the 
only surviving extra contribution to neutrino mass among those discussed is from the one 
corresponding to expression (P). Notice that the second Higgs doublet of the Zee model, 
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corresponding to H 3 here, is also assumed not to have couplings of the form LiH :i Ej. The 
condition for this LR slepton mixing contribution to be below the required m efl would be 
the same as discussed in the previous section. 

However, there is a new contribution to (m u )kk given by 

1 <K > nl(>&^f(M*,Ml ), (32) 



16tt 2 {h°f 

which is obtained from Fig. [1] with £ L replaced by h\ and h® by h® (of course with different 
couplings at the vertices.) This is in fact a consequence of the fact that the term \ H k H^H^E c k 
provides new mass mixings for the charged Higgsinos and the charged leptons. As with the 
mixings induced by the /Xj's, the new effect is see-saw suppressed; but unlike the //j's their 
magnitude may be less severely constrained. Nevertheless, the essential difference here is 
that unlike the //j terms the X^H^HfE^. term does not contribute to the mixings between 
neutrinos and the gauginos and higgsinos on tree level. 

We will assume also that the deviations of the charged-lepton mass eigenstates resulted 
from the new mixing are negligible. Bleaching the assumption actually does not cause too 
much trouble though. Its main effect is simply the modification of the numerical values of 
rn e . , s used as the physical masses get extra contributions. Here, similar to the above we are 
interested in only the minimal set of couplings {A 12 fc , A 13 ^ , Af} with a specific k. For 
expression (|3lD to give the right value to m efJi , we need 

Max(M, 2 _,M? ) (h o\ 
A 12fe Af ^ W) (?X 10-7 ^V" 1 ) > ( 33 ) 

and similarly for m er , it requires A 13 fc = (rrvf l /m^.)X 12 k ■ This condition is easy to satisfy, 
for example, when we take (h®)/(h®) = 0.1. Next, we compare the expression ( |3~T| ) with Eq. 
(|32|). For Eq. (|3lD to dominate over Eq. (|32|), it is required that 



_ . in ) / / ' 

A12A: 3> A 
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The most favorable scenario under the context is obtained by taking k = 1 where m t is just 
the m e . The above requirements are then easily satisfied. In addition, the corresponding 
requirement for subdomination of the LR slepton mixing contribution discussed above is 
then the same as Eq. fll7D , and we also have Eq. ( |I9|) from the tree-level Zee-scalar in- 
duced muon decay. All these constraints can now be easily satisfied. Hence, having such a 
super symmetric Zee model looks very feasible. 

One may argue, following the spirit of the single- VEV parametrization, that H 3 may 
be arranged to have no VEV. That would apparently kill the scenario. However, we have 
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an assumption above that there is no LiH 3 E^ term in the superpotential, which is also 
adopted in the original Zee model. Without the assumption we could then switch H 3 and H x 
around, and though the scenario is still viable it, however, becomes much more complicated 
to analyze. Furthermore, if H x was the only one with a VEV, hs~ should take over the role 
of hi in Fig. [I], and then the couplings of the hf to leptons could not be taken diagonal in 
general. Studies of these more general situations, together with more admissible terms in 
the superpotential involving H a , actually worth more attention. This is however beyond the 
scope of the present paper. 

An alternative approach is to give up identifying the right-handed slepton as the Zee 
scalar. One can introduce a vectorlike pair of Zee (singlet) superfields E z and E£ with the 
scalar component of the latter as the Zee scalar. A X^LiLjE^ term takes the role of the Ay* 
above. The F term of with nonzero ^L^H 2 and Y^L^H^E^ terms provides the mixing 
between the new Zee scalar and the h\. But the Y% coupling easily messes up the identity 
of the physical charged leptons. It is clear then this is an even more complicated situation 
than the previous one, and has to be analyzed carefully in a different framework. 

Finally, one can take the trivial supersymmetrization by taking both E z and E% as well 
as H 3 and H 4 . The restrictions on the parameter space of the relevant couplings are then 
unlikely to have any interesting feature beyond that of the Zee model itself. It is interesting, 
however, to note that the couplings needed, L^L-jE^ and H X H Z E^ do not break R parity at 
all, though the lepton number is violated. 

VI. CONCLUSIONS 

We have discussed the embedding of the Zee neutrino mass model in the framework of 
.R-parity-violating supersymmetry. It is a nontrivial supersymmetrization of the Zee model 
in the sense that one or both of the extra scalar fields of the Zee model are identified 
within the minimal supersymmetric SM spectrum. We have analyzed in detail the minimal 
scheme where a right-handed slepton plays the role of the charged Zee scalar, with no extra 
Higgs doublet introduced. The A^ and //j couplings in i?-parity-violating supersymmetry 
are identified as the lepton-number-violating couplings in the Zee model. Nevertheless, the 
scheme also introduces other types of contributions to neutrino masses. We have described 
in detail all these contributions and their general forms. 

We have addressed and answered the question on the feasibility of such a model. For 
the minimal scheme we illustrated that a set of i?-parity-violating couplings given as 
{X 12 k , X 1S k , A 4 *:} with a specific k completes the model. The various contributions to 
neutrino masses are discussed and the conditions for maintaining the zeroth order Zee mass 
texture and for fitting the experimental neutrino oscillation data are derived. The case with 
k = 3 has been shown to be marginally feasible, though the constraints are very stringent. 
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Here, the right-handed stau is the Zee scalar. The analysis also illustrates the interesting 
interplay between different couplings, and the relative strength of each type of contributions. 
In particular, we have discussed two contributions both involving the bilinear .R-parity vi- 
olating couplings fi k , and the trilinear ones. Both contributions have not been discussed 
before and the Zee mechanism does correspond to one of them and the another one involves 
a new source of LR slepton mixing. 

We have also discussed alternatives to the minimal scheme. Among the alternatives, we 
have considered the interesting case with an extra pair of vectorlike Higgs superfields. A 
coupling for the superpotential term HfHfE%. replaces fi k of the previous case. Constraints 
on this extended type of models are much weaker, and so having a phenomenologically viable 
model of this type would not be a problem. The best scenario in this case is for k — 1, 
namely, taking the selectron as the Zee scalar. More detailed studies of such models worth 
a serious effort. We hope to report on that in the future. 
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FIG. 1. The two Feynman diagrams for the Zee mechanism in the i?-parity violating SUSY 
framework. The original charged singlet boson h% ee of the Zee model is shown in parentheses. 
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FIG. 2. The one- loop diagrams for neutrino mass generated by the usual term m t (A E —fi tan 0) 
in LR slepton mixing. 
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FIG. 3. The one-loop diagram for neutrino mass generated by the term (h®) fJ- g ^ g hk^L h ^* Rk in 
the LR slepton mixing. 
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